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Abstract 

Data of the Crystal Barrel and L3 collaborations clarified essentially the situation 
in the 2 ++ sector in the mass region up to 2400 MeV, demonstrating the linearity of 
(n, M 2 ) trajectories, where n is the radial quantum number of a quark-antiquark state 
with mass M. We discuss these data and show that there exists a superfluous state for the 
(n, M 2 ) trajectories: a broad resonance /2(2000). We pay special attention to the reactions 
pp — > 7T7T, rjrj, rjrj' in the mass region 1990-2400 MeV where, together with /2(2000), four 
relatively narrow resonances are seen: /2(1920), /2(2020), /2(2240), /2(2300). We analyse 
the branching ratios of all these resonances and show that only the decay couplings of the 
broad state /2(2000) — ► 7r°7r°, rjrj, rjrj 1 satisfy relations inherent in the glueball decay. 

PACS numbers: 14.40-n, 12.38-t, 12.39-MK 



1 Introduction 

A broad isoscalar-tensor resonance in the region of 2000 MeV is seen in various reactions [1]. 
Recent measurements give: 

M = 2010 ± 25 MeV, T = 495 ± 35 MeV in pp -> 7r°7r°, 7777, rjrj' [2], 
M = 1980 ± 20 MeV, T = 520 ± 50 MeV in pp -> ppirmrir [3] , 
M = 2050 ± 30 MeV, T = 570 ± 70 MeV in ir'p -> <p<pn [4]; 
following them, we denote the broad resonance as / 2 (2000). 

The large width of / 2 (2000) arouses the suspicion that this state is a tensor glueball. Such 
an opinion was expressed lately in different publications. 
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In [5], Chapter 5.4, it is said that the very broad isoscalar 2 ++ state observed in the region 
~ 2000 MeV with a width of the order of 400 - 500 MeV [3] could well be the trace of a tensor 
glueball lying on the Pomeron trajectory. 

Another argument comes from the analysis of the mass shifts of the qq tensor mesons ([6], 
Section 12). It is stated here that the mass shift between / 2 (1560) and a 2 (1700) can not be 
explained by the mixing of non-strange and strange components in the isoscalar sector: in 
such a mixing the average mass squared does not change and we should find / 2 (1750) at a 
much higher mass. Instead, we observe a shift down in masses of both isoscalar states. Such 
a phenomenon can be an indication for the presence of a tensor glueball in the mass region 
1800-2000 MeV. 

In [4], the following argument is presented: a significant violation of the OZI-rule in the pro- 
duction of tensor mesons with dominant ss components (reactions ir~p — > / 2 (2120)n, / 2 (2340)n, 
/ 2 (2410)n — » 00n [7]) is due to the presence of a broad glueball state / 2 (2000) in this region, 
resulting in a noticeable admixture of the glueball component in / 2 (2120), / 2 (2340), / 2 (2410). 

The possibility that the broad resonance / 2 (2000) could be a glueball is discussed also in 
[8], Section 10. Here, however, a problem in the identification of / 2 (2000) as the tensor glueball 
is stressed. As it is written in Subsection 10.6, of [8] the prediction for the branching fraction 
of the 2 + glueball is large if the width is taken to be the 500 MeV fitted to / 2 (1950). Observed 
decays to oo and / 2 (1270)cr account for (10 ± 0.7 ± 3.6) ■ 10~ 4 of J radiative decays and for 
a further (7 ± 1 ± 2) • 10~ 4 in K*K* decays. If one assumes flavour-blindness for vector-vector 
final states, the vector-vector contribution increases to (16 ± 2 ± 4.5) • 10~ 4 . The total 2.6 • 10~ 3 
is still less by a factor 9 than predicted for a glueball; in [8] this is considered as a problem in 
identifying / 2 (1950) with the 2+ glueball. 

In [9] it was emphasised that the / 2 (2000) being superfluous for qq systematics can be 
considered as the lowest tensor glueball. A recent re-analysis of the 00 spectra [4] in the 
reaction ir~p — > 00n [7], the study of the processes 77 — > 7r + 7r~7r° [10], 77 — > K$Ks [6] and the 
analysis of the pp annihilation in flight pp — > titt, rjr], i]i] f [2] clarified essentially the status of the 
( J PC = 2 ++ )-mesons. This allows us to place the / 2 mesons reliably on the (n, M 2 ) -trajectories 
[9], where n is the radial quantum number of the gg-state. In the present review we discuss the 
data for 77 — > 7r + 7r~7r° [10], 77 — > KsKs [6] and pp — > nn, rjr), rjrj' [2] in Section 2. 

In [11] (see also [5, 12]), the known gg-mesons consisting of light quarks (g = u,d,s) are 
placed on the (n, M 2 ) trajectories. Trajectories for mesons with various quantum numbers turn 
out to be linear with a good accuracy. In Section 3 we give a systematisation of tensor mesons, 
f'2 and a 2 , on the (n, M 2 ) planes. 

The quark states with (/ = 0, J PC = 2 ++ ) are determined by two flavour components 
nn = (uu + dd)/\^2 and ss for which two states 2S+1 Lj = 3 P 2 , 3 F 2 are possible. Consequently, 
we have four trajectories on the (n, M 2 ) plane. Generally speaking, the / 2 -states are mixtures 
of both the flavour components and the L — 1, 3 waves. The real situation is, however, such 
that the lowest trajectory [/ 2 (1275), / 2 (1580), / 2 (1920), / 2 (2240)] consists of mesons with dom- 
inant 3 P 2 nn components, while the trajectory [/ 2 (1525), / 2 (1755), / 2 (2120), / 2 (2410)] contains 
mesons with predominantly 3 P 2 ss components. The F-trajectories are presently represented 
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by three resonances [/2(2020), /2(2300)] and [/2(2340)] with the corresponding dominant 3 F 2 nn 
and 3 F 2 ss states. In [9], it is shown that the broad resonance / 2 (2000) is not part of those states 
placed on the (n, M 2 ) trajectories. In the region of 2000 MeV three nn-dominant resonances, 
/2(1920), / 2 (2000) and / 2 (2020), are seen, while on the (n, M 2 ) -trajectories there are only two 
vacant places. This means that one state is obviously superfluous from the point of view of the 
gg-systematics, i.e. it has to be considered as exotics. 

There exist various arguments in favour of the assumption that / 2 (2000) is generated by a 
glueball. Still, it can not be a pure gluonic / 2 (2000) state: it follows from the 1/N expansion 
rules [33, 34] that the gluonic state (qq) mixes with quarkonium systems (gg) without suppres- 
sion. The problem of the mixing of (gg) and (qq) systems is discussed in Section 4, where we 
present also the relations between decay constants of a glueball into two pseudoscalar mesons 
glueball — > PP and into two vector mesons glueball — ■> VV. 

In Section 5 we demonstrate that just / 2 (2000) is the glueball. 

In [13] it was pointed out that an exotic state has to be broad. Indeed, if an exotic resonance 
occurs among the standard gg-states, they overlap, and their mixing becomes possible due 
to large distance transition: the resonance (1) — > real mesons — > resonance (2) . Owing to 
these transitions, an exotic meson accumulates the widths of its neighbouring resonances. The 
phenomenon of the accumulation of widths was studied in the scalar sector near 1500 MeV 
[14, 15]. In [15], a model of mixing of the gluonium gg with the neighbouring quarkonium 
states was considered. It was demonstrated that, as a result of mixing, it is precisely the 
gluonium state which transforms into a broad resonance. The reason is that gg states being 
orthogonal to each other mix weakly, while the gluonium mixes with neighbouring gg states 
without suppression. Therefore, the gluonium "dives" more rapidly into the complex M-plane. 
The mixing of states is always accompanied by a repulsion of the corresponding poles: when 
poles are in the complex M-plane at approximately the same Re M, this repulsion results in 
"sinking" one of them into the region of large and negative Im M and "pulling" others to the 
real M-axis. 

Hence, the large width of / 2 (2000) can indicate that this state is an exotic one. Strictly 
speaking, this fact is not sufficient to prove that / 2 (2000) is a glueball. At the moment a variety 
of versions for exotic mesons is discussed; these are qqg hybrids as well as multiquark states 
(see, e.g. [5, 8, 12] an references therein). Thus, in order to fix / 2 (2000) as a glueball, it is 
of great importance to investigate the decay couplings and prove that they satisfy relations 
characterising the glueball. The coupling constants for the transitions 

/ 2 (1920), / 2 (2000), / 2 (2020), / 2 (2240), / 2 (2300) - tttt, nn, nn' 

are separated in [16, 17] on the basis of a partial wave analysis [2] carried out earlier. The 
coupling constants obtained in [16, 17] indicate that only the decays / 2 (2000) — > n n ,nn,nn' 
obey the relations corresponding to a glueball, while the decay constants for other resonances 
do not fulfil such conditions. Note that the glueball decay couplings are close to those for the 
SU (3)-Q&vour singlet, but, because of the flavour symmetry violation caused by the strange 
quark, do not coincide with them exactly. 

Let us remind that there are two more arguments in favour of the glueball nature of / 2 (2000): 
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(i) the Pomeron trajectory, determined on the basis of data on high-energy hadron decays (see, 
e.g., [18, 19, 20]), indicates that a tensor glueball has to have a mass of the order of 1.7—2.3 GeV; 

(ii) lattice calculations [21] lead to a similar value, M 2 +- g i ue b a u ~ 2.3 — 2.5 GeV. 

We have one, sufficiently general, argument against the interpretation of / 2 (2000) as an ex- 
otic qqg or qqqq state: the absence of any serious facts confirming their existence. Indeed, if such 
states existed, we would see a large number of them in the mass region above 1500 MeV. More- 
over, we could observe not only exotic states; the number of resonances with "normal" quantum 
numbers would also be seriously increased. However, the systematics of quarks on the (n, M 2 )- 
plane does not reveal such an increase: almost all observed resonances can be interpreted as 
qq states (see [5], Chapter 5). Apparently, Nature does not like coloured multiparticle objects. 
The same conclusion follows from the systematisation of baryons: experimental data give a 
much smaller amount of excited states, than calculations in the framework of a three-quark 
model [22] do. One gets the impression that excited baryons are rather quark-diquark systems 
(see discussions in [5, 23]. 

Owing to the 1/N expansion rules, the gluonium component is relatively small in the quark 
state fa: its probability is suppressed as 1/N C . In Section 6 we determine the mixing angle of 
the nn and ss components in the quark / 2 -mesons, making use of the relations between the 
decay constants fa — > 7777,7777,7777'. Also, we estimate the possible changes in the mixing angle 
as a consequence of a gluonium component in the quark state fa. 

In the Conclusion, we discuss the situation in the glueball sector. 

Up to now, two glueball states, the scalar meson /o(1200 — 1600) an the tensor ^(2000) are 
observed. Both states are broad ones, and the coupling constants corresponding to their decays 
into pseudoscalar mesons (channels glueball — > PP) satisfy just the relations characterising the 
glueball. The next states which are of interest are radial excitations of the scalar and tensor 
gluonia, and the pseudoscalar glueball. Taking the Pomeron trajectories on the (J, M 2 )-plane 
as a basis, we predict the masses of excited scalar and tensor glueball states. 

2 Analysis of the data for tensor mesons 

We demonstrate here the results obtained from the data analysis used in performing the sys- 
tematisation of the fa resonances and extracting the decay couplings fa — > 7777,7777,7777'. 

2.1 L3 data on the 77 — > 77+77 77° reaction 

In this reaction the 77 channel couples only to states with C=+l parity; 377 has a negative 
G-parity. For a qq system one has G = Ce t7Tl , so the 1 = 1 quark-antiquark states are produced 
only in the 77 channel. Due to C-parity conservation in neutral decay modes, only /-states 
with (J PC = ++ , 2 ++ , 4 ++ ...) are produced in the 77+77" channel. In the 77+77° channel only 
isovector mesons with J PC = 1 , 3 . . . are produced. 

The 77 mass distribution is dominated by the production of the a2(1320) resonance, see 
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Figure 1: Reaction 77 — > 7r + 7r~7r°. a) 77 mass spectrum, b,c) 7r + 7r~ and -k ± tt mass spectra, 
d,e) the angular distributions of charged and neutral pions in the c.m.s. of the reaction, f) the 
angular distribution between charged and neutral pion in the c.m.s. of two charged pions. 



Fig. la. One can see a prominent structure in the mass region 1.6-1.8 GeV as well as a possible 
contribution at the a-i (1320) signal. 

The 7r + 7r~ mass distribution is shown in Fig. lb. There are no clear signals in the data 
coming from well known narrow scalar-isoscalar states / (980) and / (1500). Indeed, the partial 
wave analysis shows very small contributions of these mesons; such decay modes were omitted 
in the final fit. 

A signal coming from / 2 (1275)7r° is observed at high 77 energies; this is important to 
describe the two-pion mass spectrum and angular distributions. 

The 7T7T — > 7T7T S-wave amplitude has a broad component which covers the mass region from 
the mi threshold up to 2 GeV. Such a component is introduced in the present analysis and is 
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Figure 2: Reaction 77 — > 7r + 7r~7r°. The contribution of different resonances listed in Table 1 to 
the cross section: full curves - 2 ++ states, dashed curves - CT + states and dotted curves - the 
contribution of 2~ + states. 

parametrized in two different ways. 

The first parametrisation is taken from [24]. It was introduced to describe the CERN-Munich 
data on (n~p — > n + 7i^n) [25] and the Crystal Barrel data on proton-antiproton annihilation 
into 3tc° and 2r]7r° channels simultaneously. To simulate a possible s-dependence of the vertex 
(which can be important for this very broad state), we use the method suggested in [24]. 

The second parametrisation was used in [26]. It covers the mass region up to 1.9 GeV 
and describes, in the framework of the P-vector/K- matrix approach, a much larger number of 
two- and three-body reactions. To avoid an over-parametrisation of the fit, we vary only the 
production couplings of the two lowest K-matrix poles. 

The main signal in the n^n mass spectrum is due to the production of p(770). There is 
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very little structure in the region higher than 1 GeV (see Fig. lc) and the signal is almost zero 
at masses above 1.5 GeV. Because of this, neither p 3 (1690), nor p(1770) has to be introduced in 
the analysis. A contribution from p(1450) is found to be useful to describe the data: however, 
this state is quite broad and possibly simulates a non-resonant two pion production in this 
channel. 

In the 7r + 7r~7r° spectrum one can see a strong signal coming from a2(1320); the characteristics 
of this resonance were defined with high precision by the VES collaboration [27]. It is not 
surprising that the 77 — > 3n data are dominated by the production of the a 2 (1320) state, since 
this resonance has the highest spin in the mass region below 1.6 GeV (the 77 cross section is 
proportional to (2 J +1)) and it is a ground qq state with the radial quantum number n — 1. 
Indeed, the 77 — > resonance production amplitude is a convolution of the photon and the 
quark-antiquark resonance wave functions [28, 29, 30]. This provides a suppression of nearly 
an order for the production of the radially excited states (n > 2) [29]. Nevertheless, there is 
a manifest contribution of the higher tensor state. While a 2 (1320) decays practically only into 
the p(770)-7T channel, the second tensor state decays almost equally to p(770)7r and / 2 (1275)7r. 
This fact allows us to identify this state with a good accuracy. 

The solution reveals quite a large contribution coming from the 0~ + partial wave decaying 
into JqIT. There is, however, a problem in distinguishing it from the experimental background: 
such a decay, giving S-wave amplitudes in all decay channels, provides very smooth structures in 
mass distributions and, moreover, these amplitudes do not interfere with the tensor amplitude. 

The contribution of the 2~ + state is found to be quite small when fitted to the 7r 2 (1670) 
state (see Table 1). If it is fitted with free Breit-Wigner parameters, it is always optimised at 
higher masses (~ 1870 MeV). 

The ir + -K~n° spectrum and the contributions of resonances with different J PC in the final 
solution are shown in Fig. 2. Masses, total widths, the T 77 partial width and the branching 
ratio into 3ir are listed in Table 1 for the considered resonances. 



Table 1: Masses, total widths and the partial width T^^ttk for the resonances observed in the 
reaction 77 — > 7r + 7r~7r°. 



Resonance 


M (MeV) 


T(MeV) 


r 77 Br(37r)(keV) 


a 2 (1320) 
a 2 (1730) 
tt(1300) [1] 

2-+ 

tt 2 (1670)[1] 


1302 ± 3 ± 6 
1725 ± 25 ± 10 
1350 ± 40 
1870 ± 60 
1670 


118 ±6 ± 10 
340 ± 40 
320 ± 50 
325 ± 40 
260 


0.65 ±0.05 

o.34i8:JS 

< 0.8 

0.15 ±0.03 

< 0.1 
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2.2 L3 data on the reaction 77 — > KsKs [6] 

Important data for the identification of the tensor mesons were obtained by the L3 collaboration 
on the reaction 77 — > K$Ks- Only states with even spin J and positive parities P = C = 
+ contribute to two neutral pseudoscalar particles, what reduces the possible partial waves 
drastically. The (2 J + 1) factor in the cross section favours the dominant production of the 
tensor states. The 4 ++ states are produced at high energies (M > 1.9 GeV). 

The meson states consisting of light quarks u, d s form meson nonets: three isospin 1, four 
isospin 1/2 and two isoscalar states. The isoscalar states can be a mixture of nn = (uu+dd)/y/2 
and ss components. The decay of (I = 0) and (I = 1, 1 3 = 0)-states into two kaons is defined by 
the production of a new ss pair (an s-quark exchange process) and has the following structure for 
different isospins: (uu+dc£)/y/2 -> K + K-+K°K° for 1=0, and (uu-dd)/V^ -> K + K--K°K° 
for 1=1. As a result, a strong destructive interference occurs between the /2(1275) and ci2(1320) 
mesons which is very sensitive to the mixing angle of the isoscalar state. The flavour content 
of isoscalar-scalar resonances belonging to the same nonet can be written in the form 



Although the production of states with dominant ss components is suppressed by the smaller 
77 coupling, these states usually have a larger branching ratio for the decay into the K$Ks 
channel, and can contribute appreciably into the total cross section. 

There is no doubt about the nature of tensor resonances below 1600 MeV. The partial wave 
analysis showed that tensor resonances are produced dominantly in the 5 S , 2 77 state, which was 
predicted by model calculations [29]. The model gives directly the ratio between 5 S 2 and 1 _D 2 
waves and the 77 couplings. These values can be introduced in the analysis without changing 
the quality of the description. The data with the lowest tensor states are shown in Fig. 3a. 

We have found that the ++ partial wave can be fitted well in the framework of the P- 
vector/K- matrix approach. The form of the ++ contribution in 77 — > KK follows closely the 
7T7T — > KK cross section, which is not surprising: the production coupling to the ss system is 
strongly suppressed in both reactions. The values of the coupling constants of the reactions 
77 — > fo,a agree well with those given in the calculations [29]: in the same way as in the 
tensor sector, the calculated 77 couplings can be used directly, not damaging the quality of the 
description. 

A contribution of 4 ++ states is observed in high energy angular distributions. Some broad 
and some rather narrow components are seen in this wave. The broad state can be associated 
with nn and the narrow one with a 4 ++ ss state. The description of the data with tensor, 
scalar and 4 ++ states is shown in Fig. 3b. 

There is a clear resonance structure in the 1750 MeV region. The angular distribution in this 
region follows very closely the (1 — cos 2 ©) 2 shape which corresponds to the 5 S , 2 77 production 
of the tensor meson. However, the acceptance decreases rapidly in the forward and backward 
directions providing a very similar dependence. Due to such a behaviour and to the not too 
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high statistics, the partial wave analysis produces almost the same angular distribution for 
a tensor state and for a scalar state. Still, the fit using a scalar state fails to reproduce the 
structure in the 1700-1800 MeV mass region. The description of the data with the best x 2 is 
shown in Fig. 3c. The mass was optimised to 1805 ± 30 MeV and the width to 260 ± 30 MeV. 
With such a mass and width, the /o state can describe the slope in mass distribution above 
1800 MeV. If the mass and width of the scalar state are fixed at the BES result [31] M = 1740 
MeV, we obtain a description shown in Fig. 3d. 

The main problem in a fit with a /o state is that there is no way to reproduce the dip 
in the 1700 MeV region and the slope above 1800 MeV using any (even very sophisticated) 
parametrisation of the resonance width. A f state can interfere only with the 1 D 2 component 
of a 2 ++ state, and this partial wave is very small in the data. Consequently, /o and / 2 
contributions practically do not interfere and it is impossible to create a dip in the data. 
Contrary to this, a tensor resonance interfering with the tails of other tensor states naturally 
produces a dip and a good description of this mass region. 

All isoscalar and isotensor states can contribute to the 77 — > K S K S cross section; this situ- 
ation offers a very good possibility to study the reaction on the basis of the nonet classification. 
With SU(3) relations imposed, the only parameters to fit the data for the first three states are 
masses, widths, the mixing angle and SU (3) violation factors. We found all masses and widths 
for the members of the first tensor nonet to be in a very good agreement with PDG. 

To describe the second nonet, we fixed parameters for / 2 (1560) from the Crystal Barrel 
results and for a 2 (1700) from the latest L3 analysis of the 7r + 7r~7r° channel [10]. At a given 
mixing angle the nonet coupling was calculated to reproduce the nn width of 20-25 MeV for the 
/ 2 (1560) state. We found a very good description of the data with SU(3) relations imposed, 
similar to the fit described in the previous section. The masses, widths, radii, KK couplings, 
mixing angles and partial widths of the states are given in Table 2; the description of the data 
is shown in Fig. 4. 

The / 2 (1750) state has a mass 1755 ± 10 MeV and a total width 67± 12 MeV. Its KK width 
is 23 ± 7 MeV; the rest of the width is likely to be defined by the K*K channel. This resonance 
destructively interferes with the tail of the / 2 (1525) state, creating a dip in the mass region 
1700 MeV. When the sign of the real part of the / 2 (1750) amplitude changes, this interference 
becomes positive, producing a clear peak in the data. 

A serious problem appears in the description of the data in the framework of the nonet 
approach, if the peak at 1750 MeV is assumed to be owing to a scalar state. If this state is a 
nonet partner of one of the known states, e.g. /o(1370) or / (1500), then the calculated signal 
is too weak to fit the data. If the KK coupling of this scalar state is fitted freely, we find that 
it must be about four times larger than the total width of the resonance. This is due to the 
2J+ 1 suppression factor and to the absence of a positive interference with the tail of / 2 (1525) 
which boosts the peak in the case of a tensor state. These are problems additional to those 
connected with the description of the dip in the 1700 MeV region. 
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Figure 3: Reaction 77 — > KsKs'- a) the description of the data with the first tensor nonet, 
b) same as (a) plus scalar states and a 4 ++ state, c) same as (b) but the scalar states with 
free parameters , d) same as (b) but / (1710) state with parameters fixed from the latest BES 
results. 
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Figure 5: Angular distributions in the reactions pp — > im, r)T), r)T)' and their fit within resonances 
ofEq.(2). 
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Figure 6: Differential cross sections in the reaction pp — > 7r + 7r [42] at proton momenta 360- 1 300 
MeV and their fit within resonances of Eq. (2). 
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Figure 7: Differential cross sections in the reaction pp — > n + n [42] at proton momenta 1350- 
2230 MeV and their fit within resonances of Eq.(2). 
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Table 2: Parameters of resonances observed in the reaction 77 — > K$Ks- The values marked 
with stars are fixed by using other data. 





First nonet 


Second nonet 




a 2 (1320) 


/ 2 (1270) 


£(1525) 


a 2 (1730) 


/ 2 (1560) 


/ 2 (1750) 


Mass (MeV) 


1304 ± 10 


1277 ± 6 


1523 ± 5 


1730* 


1570* 


1755 ± 10 


Width (MeV) 


120 ± 15 


195 ± 15 


104 ± 10 


340* 


160* 


67 ± 12 


9 L (GeV) 


0.8 ±0.1 


0.9 ±0.1 


1.05 ±0.1 


0.38 ±0.05 


up (deg) 


-1 ±3 


1U -10 




Figure 8: Polarisation in pp — > n + n [42] and its fit within resonances of Eq. (2). 
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Figure 9: Cross sections and Argand-plots for 3 P 2 and 3 F 2 waves in the reaction pp — > 
7r°7r°, 7777, 7777'. The upper row refers to pp — > 7r°7r°: we demonstrate the cross sections for 
3 P 2 and 3 F 2 waves (dashed and dotted lines, respectively) and the total (J = 2) cross section 
(solid line) as well as Argand-plots for the 3 P 2 and 3 F 2 wave amplitudes at invariant masses 
M = 1.962, 2.050, 2.100, 2.150, 2.200, 2.260, 2.304, 2.360, 2.410 GeV. The figures in the second 
and third rows refer to the reactions pp — > r\r\ and pp — > r/r]'. 
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2.3 Data for proton-antiproton annihilation in flight pp — > 7T7t, 7777, 7777' 

The pp annihilation in flight gives information about resonances with M > 1900 MeV. High 
statistical data taken at antiproton momenta 600, 900, 1150, 1200, 1350, 1525, 1640, 1800 and 
1940 MeV/c were used for the analysis pp — > 7r°7r°, rji], rjr]' [2]. The combined analysis was 
performed together with pp — > 7r + 7r~ data obtained with a polarised target [42]. Five tensor 
states are required to describe the data, / 2 (1920), / 2 (2000), } 2 (2020), / 2 (2240), / 2 (2300): 

Resonance Mass(MeV) Width(MeV) (2) 

/ 2 (1920) 1920 ±30 230 ±40 

/ 2 (2000) 2010 ±30 495 ±35 

/ 2 (2020) 2020 ±30 275 ± 35 

/ 2 (2240) 2240 ±40 245 ± 45 

/ 2 (2300) 2300 ±35 290 ± 50 . 

The description of the data is illustrated by Figs. 5,6,7 and 8. In Fig. 9 we show the cross 
sections for pp — > tt°tt , rji], rji]' in 3 P 2 pp and 3 -F 2 pp waves (dashed and dotted lines) and the total 
(J = 2) cross section (solid line) as well as Argand-plots for the 3 P 2 and 3 F 2 wave amplitudes 
at invariant masses M = 1.962, 2.050, 2.100, 2.150, 2.200, 2.260, 2.304, 2.360, 2.410 GeV. 

The pp — > 7r°7r°, 7777, 777/ amplitudes provide the following ratios for the / 2 resonance cou- 
plings g^o : g m : g m r. 

/ 2 (1920) 1 : 0.56 ±0.08 : 0.41 ±0.07 (3) 

/ 2 (2000) 1 : 0.82 ±0.09 : 0.37 ±0.22 

/ 2 (2020) 1 : 0.70 ±0.08 : 0.54 ±0.18 

/ 2 (2240) 1 : 0.66 ±0.09 : 0.40 ±0.14 

/ 2 (2300) 1 : 0.59 ±0.09 : 0.56 ±0.17. 

These coupling ratios allow one to estimate the quarkonium-gluonium c content of the / 2 - 
resonances. 



3 Systematisation of tensor mesons on the (n, M 2 ) tra- 
jectories 

In [11] (see also [5, 12]), the known gg-mesons consisting of light quarks (q = u, d, s) were put 
on the (n, M 2 ) trajectories, where n is the radial quantum number of the qq system with mass 
M. The trajectories for mesons with various quantum numbers turn out to be linear with a 
good accuracy: 

M 2 = M 2 ±(n-l)/i 2 (4) 

where /z 2 = 1.2 ± 0.1 GeV 2 is a universal slope, and M is the mass of the lowest state with 
n — 1. 



17 



CM 

QJ 
<3 



1 a ) / 






2470 ./ 


: few ^ 300 


X F2240 




2720. 




- f2020 / 






1 1920 


1755y( / 




; i525 ./ /* 1580 




</l285 




: / 2 (0 + 2 ++ ) 

. . . 1 . . . 1 . . . 


\l 2 =1.15GeV 2 ' 
. . . 1 . . . 1 . . . 



6 
n 



"M 6 

(3 

CM- 5 



b) 




2255 



1320 



a 7 (l'2 ++ ) 



\l 2 =1.15GeV 2 



4 



6 
n 



Figure 10: The f 2 and a 2 trajectories on the (n, M 2 ) plane; n is the radial quantum number of 
the qq state. The numbers stand for the experimentally observed / 2 and a 2 masses (M). 



In Fig. 10a we demonstrate the present status of the (n, M 2 ) trajectories for the f 2 mesons 
(i.e. we use the results given by [2, 4, 6]). To avoid confusion, we list here the experimentally 
observable masses. First, it concerns the resonances seen in the 00 spectrum [7]. In [4] the 00 
spectra were re-analysed, taking into account the existence of the broad / 2 (2000) resonance. As 
a result, the masses of three relatively narrow resonances are shifted compared to those given 
in the compilation PDG [1]: 

/ 2 (2010)| PDG -> / 2 (2120) [4], / 2 (2300)| PDG -> / 2 (2340) [4], / 2 (2340)| PDG -> / 2 (2410) [4]. 

The trajectory for the a 2 -mesons, Fig. 10b, is drawn on the basis of the recent data [51]. 

The quark states with (7 = 0, J PC = 2 ++ ) are determined by two flavour components nn and 
ss for which two states 2S+1 Lj = 3 P 2 , 3 F 2 are possible. Consequently, we have four trajectories 
on the (n, M 2 ) plane. Generally speaking, the / 2 -states are mixtures of both the flavour com- 
ponents and the L = 1, 3 waves. The real situation is, however, such that the lowest trajectory 
[/ 2 (1275), / 2 (1580), / 2 (1920), / 2 (2240)] consists of mesons with dominant 3 P 2 nn components 
(we denote nn = (uu + dd)jj%, while the trajectory [/ 2 (1525), / 2 (1755), / 2 (2120), / 2 (2410)] 
contains mesons with predominantly 3 P 2 ss components, and the F-trajectories are represented 
by three resonances [/ 2 (2020),/ 2 (2300)] and [/ 2 (2340)] with the corresponding dominant 3 F 2 nn 
and 3 F 2 ss states. Following [9], we can state that the broad resonance / 2 (2000) is not part of 
those states placed on the (n, M 2 ) trajectories. In the region of 2000 MeV three nn-dominant 
resonances, / 2 (1920), / 2 (2000) and / 2 (2020), were seen, while on the (n, M 2 )-trajectories there 
are only two vacant places. This means that one state is obviously "superfluous" from the 
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point of view of the gg-systematics, i.e. it has to be considered as exotics. The large value of 
the width of the / 2 (2000) strengthen the suspicion that, indeed, this state is an exotic one. 

4 Quarkonium and qluonium states: mixing and decay 

On the basis of the 1/iV-expansion rules, we estimate here effects of mixing of quarkonium and 
qluonium states. Then, making use of the rules of quark combinatorics, we give the relations 
for decay constants of these states. 

4.1 Mixing of qq and gg states 

The rules of the 1/iV-expansion [33, 34], where N = N c = Nf are numbers of colours and light 
flavours, provide a possibility to estimate the mixing of the gluonium (gg) with the neighbouring 
quarkonium states (qq). 

The admixture of the gg component in a gg-meson is small, of the order of 1/N C : 

h(qq — meson) = gg cos a + gg sin a (5) 
sin 2 a ~ 1 /N c . 

The quarkonium component in the glueball should be larger, it is of the order of Nf/N c : 

/ 2 (glueball) = gg cos 7 + (qq) g i ue baii sin 7, (6) 
sin 2 7 ~ N f /N c , 

where (qq) g i U ebaii is a mixture of nn = (uu + dd)/\/2 and ss components: 

(qq)glueball = nn COS (f glueball + SS sin <f g l uebaU , (7) 

with sirup giuebaii = ^J\/(2 + A). If the flavour SU(3) symmetry were satisfied, the quarkonium 
component (qq)gi ueba ii would be a flavour singlet, <f g i ue baU -> fsingiet - 37°._ In reality, the 
probability of strange quark production in a gluon field is suppressed: uu : dd : ss = 1 : 1 : A, 
where A ~ 0.5 — 0.85. Hence, (qq) g i U ebaii differs slightly from the flavour singlet, it is determined 
by the parameter A as follows [43] : 

(qq) giuebaii = (uu + dd + V\ ss) I V2 + A . (8) 

The suppression parameter A was estimated both in multiple hadron production processes [44] , 
and in hadronic decay processes [45, 26]. In hadronic decays of mesons with different J PC the 
value of A can be, in principle, different. Still, the analyses of the decays of the 2 ++ -states 
[45] and ++ -states [26] show that the suppression parameters are of the same order, 0.5-0.85, 
leading to 

Vglueball - 26° - 33°. (9) 
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Figure 11: Examles of the diagrams which determine gluonium (gg) decay. 
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Let us explain now Eqs. (5)- (8) in detail. 

First, let us evaluate the transition couplings using the rules of 1/iV-expansion; this eval- 
uation will be done for the decay transitions gluonium — > two qq-mesons and quarkonium-^ 
two qq-mesons. For this purpose, we consider the gluon loop diagram which corresponds to 
the two-gluon self-energy part: gluonium — > two gluons — * gluonium (see Fig. 11a). This loop 
diagram B (gluonium — > gg — > gluonium) is of the order of unity, provided the gluonium is a 
two-gluon composite system: B (gluonium — > gg — > gluonium) ~ 9gi U onium-*gg-^c ~ 1> where 

9 gluonium^ 

gg is a coupling constant for the transition of a gluonium to two gluons. Therefore, 

g gluonium^ gg ~ l/^c • (10) 

The coupling constant for the gluonium — > gg transition is determined by the diagrams of Fig. 
lib type. A similar evaluation gives: 

<^_^ ~ ^ium^p ~ W ■ (11) 

Here g is the quark-gluon coupling constant, which is of the order of 1/ \/W c [33]. The coupling 
constant for the gluonium — > two qq-mesons transition in the leading 1/N C terms is governed 
by diagrams of Fig. 11c type: 

9 gluonium— >two mesons ~ 9 gluonium— ^qq 9meson—>qqNc ~ l/^c • (12) 

In (12) the following evaluation of the coupling for transition qq — meson — > qq has been used: 

9meson^qq ~ 

(13) 

which follows from the fact that the loop diagram of the gg-meson propagator (see Fig. 12a) is 
of the order of unity: B(qq — meson — > qq — > meson) ~ 9 me son^qqNc ~ 1 • 

The diagram of the type of Fig. lid governs the couplings for the transition gluonium — > 
£u>o qq-mesons in the next-to- leading terms of the l/iV c -expansion: 

9 gluonium^two mesons ~ 9gluonium—*gg 9meson—>gg c 

-1/7V 2 , (14) 
where the coupling g meson ^ g g has been estimated following the diagram in Fig. 12b: 



imeson^gg ymeson^qq 9 n ~ l 1/JVf. (15) 



Decay couplings of gg- meson into two mesons in leading and next-to-leading terms of 1/N C 
expansion are determined by diagrams of the type of Figs. 12c and 12d, respectively. This 
gives 



a t ~ a -N„ ~1/\ N„ (16) 

t) meson^two mesons Hmeson^qq c ^/ y J 'o V- 1 -"/ 

9 meson^two mesons ~ 9meson^qq9meson^gg9 Jv c ~ 1 / iv c 

Now we can estimate the order of the value of sin 2 7 which defines the probability (qq) g i U ebaii, 
see Eq. (6). This probability is determined by the self-energy part of the gluon propagator 
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(diagram in Fig. lie) — it is of the order of Nf/N c , the factor Nf being the light flavour number 
in the quark loop. Let us emphasise that the diagram in Fig. lie stands for only one of the 
contributions of that type; indeed, contributions of the same order are also given by diagrams 
with all possible (but planar) gluon exchanges in the quark loop. 

One can also evaluate sin 2 7 using the transition amplitude gluonium — > quarkonium (see 
Fig. 12e), which is of the order of l/yfN~ c . The value sin 2 7 is determined by the transition 
amplitude squared, summed over the flavours of all quarkonia, thus resulting in Eq. (6). 

The probability of the gluonium component in the quarkonium, sin 2 a, is of the order of the 
diagram in Fig. 12f, ~ 1/N C , giving us the estimate (5). Here, as in the self-energy gluonium 
block, planar-type gluon exchanges are possible. Because of this, in the intermediate qq state 
all the interactions are taken into account. 

The diagram in Fig. lie defines also the flavour content of (qq) g i ue bau — we see that the 
gluon field produces light quark pairs with probabilities uu : dd : ss = 1 : 1 : A, so (qq) g i U ebaii is 
determined by Eq. (8) not being a flavour singlet. 

4.2 Quark combinatorial relations for decay constants 

The rules of quark combinatorics lead to relations between decay couplings for mesons which 
belong to the same SU(3) nonet. The violation of the flavour symmetry is taken into account 
by introducing a suppression parameter A for the production of the strange quarks by gluons. 

In the leading terms of the 1/N expansion, the main contribution to the decay coupling 
constant comes from planar diagrams. Examples of the production of new gg-pairs by inter- 
mediate gluons are shown in Figs. 13a and 12b. When an isoscalar gg-meson disintegrates, 
the coupling constants can be determined up to a common factor, by two characteristics of a 
meson. The first is the quark content of the gg-meson, gg = nn cos tp + ss sin tp, the second is 
the parameter A. Experimental data provide the following values for this parameter: A ~ 0.5 
[44] in central hadron production in high-energy hadron-hadron collisions, A = 0.8 ± 0.2 [45] 
for the decays of tensor mesons and A = 0.5 — 0.9 [26] for the decays of ++ mesons. 

Let us consider in more detail the production of two pseudoscalar mesons P1P2 by f 2 - 
quarkonium and /2-gluonium: 

f 2 (quarkonium) — > hit , KK , 7777 , 777/ , 7/7/ ( 1 7) 

ji (gluonium) — > irir , KK , 7777 , 7777' , rj'r)' . 

The coupling constants for the decay into channels (17), which in the leading terms of the 1/N 
expansion are determined by diagrams of the type shown in Fig. 13, may be presented as 

g L (qq^ Pl P 2 ) = Cf^tp, \)g$, , (18) 

where Cp P (<p, A) and Cp g lP2 (X) are wholly calculable coefficients depending on the mixing angle 
ip and parameter A; g P and G P are common factors describing the unknown dynamics of the 
processes. 
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Dealing with processes of the Fig. 13b type, one should bear in mind that they do not 
contain {qq) quarkonium components in the intermediate state but {qq) continuous spectrum only. The 
states (qq)quarkonium in this diagram would lead to processes of Fig. 13c, namely, to a diagram 
with the quarkonium decay vertex and the mixing block of gg and qq components. All these 
sub-processes are taken into account separately. 

The contributions of the diagrams of the type of Fig. lid and 12d, which give the next-to- 
leading terms, g NL {qq — > -P1-P2) and g NL (gg — > -P1-P2), may be presented in a form analogous 
to (18). The decay constant to the channel P\P 2 is a sum of both contributions: 

g\qq^P x P 2 ) + g NL (qq - P ± P 2 ), (19) 
g L (gg^P x P*) + g NL (gg -> P,P 2 ). 

The second terms are suppressed compared to the first ones by a factor N c ; the experience in 
the calculation of quark diagrams teaches us that this suppression is of the order of 1/10. 

Coupling constants for gluonium decays, g L {gg — > P1P2) and g NL {gg — > P1P2), are presented 
in Table 3 while those for quarkonium decays, g L (qq — > P1P2) and g NL (qq — > P1P2), are given 
in Table 4. 

In Table 5 we give the couplings for decays of the gluonium state into channels of the vector 
mesons: gg — > ViV 2 . 



24 




Figure 13: Examples of planar diagrams responsible for the decay of the gg-state (a) and the 
gluonium (b) into two gg-mesons (leading terms in the 1/N expansion), c) Diagram for the 
gluonium decay with a pole in the intermediate gg-state: this process is not included into the 
gluonium decay vertex. 

Table 3 

Coupling constants of the /2-gluonium decaying to two pseudoscalar mesons, in the leading 
and next-to leading terms of 1/N expansion. O is here the mixing angle for 77 — 77' mesons: 
i] = nn cos 6 — ss sin 9 and rj' = nn sin 6 + ss cos 9. 
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Channel 


Gluonium decay 
couplings in the 
leading term of 
1/N expansion. 


Gluonium decay 
couplings in the 
next-to-leading term 
of 1/N expansion. 


Iden- 
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(sinG + y^cosG) 


1 


rfrf 


G L (sin 2 G + A cos 2 G) 


2G JVi (sinG + Y^fcosG) 2 


1/2 
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Table 4 

Coupling constants of the /2-quarkonium decaying to two pseudoscalar mesons in the leading 
and next-to-leading terms of the 1/N expansion. The flavour content of the /2-quarkonium is 
determined by the mixing angle ip as follows: fi{qq) = nn cosy? + ss siny? where 

nn = {uu + dd) / y/2. 





Decay couplings of 
quarkonium 

IT) 1 f 1 1 T 1 ()' Tfi'FYYl 
111 It-dU-lllg, LjCIIII 
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Decay couplings of 
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V2g NL 


(cos — ^| sin 0) x 




\/A sin y? sin 2 @) 


(cos <p cos — sin (p sin 0) 


7)7)' 


g L sin© cos© (cosy?/\/2— 




(cos© - y / |sin©)x 




\/A siny?) 


(cos 


ip sin + sin y? cos 0) 






+ (sin© + y / |cos0)x 






(cos y? sin — sin ip cos 0)] 


7} f 7}' 


5f L (sin 2 cosy9/ v / 2+ 


V2g NL 


(sin© + y / |cos©)x 




\/A sin tp cos 2 


(cos y? cos + sin ip sin 0) 



27 



Table 5 

The constants of the tensor glueball decay into two vector mesons in the leading (planar 
diagrams) and next-to-leading (non-planar diagrams) terms of 1/iV-expansion. The mixing 
angle for u — mesons is defined as: u = nn cos (py — ss sin ip v , <p = nn sin ip v + ss cos ipy. 
Because of the small value of ipy, we keep in the Table only terms of the order of ipy. 





Constants for 


Constants for 


Identity factor 




glueball decays in 


glueball decays in 
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Channel 
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5 The broad state ^2(2000): the tensor glueball 



In the leading terms of l/iV c -expansion we have definite ratios for the glueball decay couplings. 
The next-to-leading terms in the decay couplings give corrections of the order of l/N c . Let 
us remind that, as we see in the numerical calculations of the diagrams, the 1/N C factor leads 
to a smallness of the order of 1/10, and we neglect them in the analysis of the decays ji — > 



7r 7r , rjrj, rjrj'. 



Considering a glueball state which is also a mixture of the gluonium and quarkonium com- 



(20) 



ponents, we have (p — >■ (fgiuebaii = sin 1 y A/ (2 + A) for the latter. So we can write 

g L ((qq) glueball -> PiP 2 ) = 9^99 -» P1P2) 
9 L {{qq) g iuebaii ~* P[Pi) g L (99 - P[Pi) 



Then the relations for decay couplings of the glueball in the leading terms of the 1/iV-expansion 
read: 

G L 

glueball ^glueball 



v / 2~+A 



g glueball = + Asill 2 0) 

y2 + A 



alueball = G^all ^ _ ^ ^ @ cog Q ( 2 1) 



Hence, in spite of the unknown quarkonium components in the glueball, there are definite 
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relations between the couplings of the glueball state with the channels 7r°7r°, 7777, rjr)' which can 
serve as signatures to define it. 



5.1 Ratios between coupling constants of / 2 (2000) — > 7r°7r°, rjrj, r\r\' 
as indication of a glueball nature of this state 

Eq. (21) tells us that for the glueball state the relations between the coupling constants are 
1 : (cos 2 6 + A sin 2 0) : (1 - A) cos sin 0. For (A = 0.5, = 37°) we have 1 : 0.82 : 0.24, and 
for (A = 0.85, = 37°), respectively, 1 : 0.95 : 0.07. Consequently, the relations between the 
coupling constants g^o : g m : g m i for the glueball have to be 

2 ++ glueball g^o : g m : g w = 1 : (0.82 - 0.95) : (0.24 - 0.07). (22) 

It follows from (3) that only the coupling constants of the broad / 2 (2000) resonance are inside 
the intervals: 0.82 < g vv / g^n < 0.95 and 0.24 > g m > j 'g^o^o > 0.07. Hence, it is just this 
resonance which can be considered as a candidate for a tensor glueball, while A is fixed in 
the interval 0.5 < A < 0.7. Taking into account that there is no place for ^(2000) on the 
(n, M 2 ) -trajectories (see Fig. 10), it becomes evident that indeed, this resonance is the lowest 
tensor glueball. 
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Figure 14: Position of the /2-poles in the complex-M plane: states with dominant ^P^nn- 
component (full circle), 3 F2nn-component (full triange), 3 P2<ss-component (open circle), 3 i*2SS- 
component (open triangle), glueball (open square). 



5.2 Mixing of the glueball with neighbouring gg-resonances 

The position of the /2-poles on the complex M-plane is shown in Fig. 14. We see that the 
glueball state / 2 (2000) overlaps with a large group of gg-resonances. This means that there is 
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a considerable mixing with the neighbouring resonances. The mixing can take place both at 
relatively small distances, on the quark-gluon level (processes of the type shown in Fig. 12e), 
and owing to decay processes 



f2(glueball) — > real mesons — > f2(qq — meson). 



(23) 



Processes of the type of (23) are presented in Fig. 15. 
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Figure 15: Transitions f 2(1^1) — > real mesons — * J '2(^2), responsible for the accumulation of 
widths in the case of overlapping resonances. 

The estimates which were carried out in Section 4 demonstrated that even the mixing at 
the quark-gluon level (diagrams of the types in Fig. 12e) leads to a sufficiently large admixture 
of the quark- ant iquark component in the glueball: f2(glueball) = cos^gg + sin^gg with sin 7 ~ 



sjNf/N c . A mixing due to processes (23), apparently, enhances the quark- ant iquark component. 

The main effect of the processes (23) is, however, that in the case of overlapping resonances one 
of them accumulates the widths of the neighbouring resonances. The position of the /2-poles 
in Fig. 14 makes it obvious that such a state is the tensor glueball. 

A similar situation was detected also in the sector of scalar mesons in the region 1000 — 1700 
MeV: the scalar glueball, being in the neighbourhood of gg-resonances, accumulated a relevant 
fraction of their widths and transformed into a broad / (1200 — 1600) state. Such a transfor- 
mation of a scalar glueball into a broad state was observed in [12, 26]; further investigations 
verified this observation. The possibility that a scalar (and tensor) glueball may considerably 
mix with gg-states was discussed already for quite a long time, see, e.g., [32, 35, 36]. At the 
same time there is a number of papers, e.g. [35, 36], in which the mixing due to transitions (23) 
is not taken into account. Hence, in these papers relatively narrow resonances like / (1500) 
and / (1710) are suggested as possible scalar glueballs. 

We see that both glueballs, the scalar / (1200 — 1600) and the tensor / 2 (2000) one, reveal 
themselves as broad resonances. We can suppose that this is not accidental. In [15] the 
transition of the lowest scalar glueball into a broad resonance was investigated in the framework 
of modelling decays by self-energy quark an gluon diagrams. As it was discovered, it was just the 
glueball which, appearing among the gg-states, began to mix with them actively, accumulating 
their widths. Hence, the glueball turned out to be the broadest resonance. 
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5.2.1 The mixing of two unstable states 



In the case of two resonances, the propagator of the state 1 is determined by the diagrams of 
Fig. 16a. With all these processes taken into account, the propagator of the state 1 is equal to: 



D 11 (s) = Iml-s-Bnis)- 



B 12 (s)B 2l (s) 
ml - s - B 22 (s) / 



(24) 



Here mi and m 2 are masses of the input states 1 and 2, and the loop diagrams B^s) are 
defined by the spectral integral 



Bi 



'<*)=/- 

J 7T 

4m 2 



d(s')g t (s') gj (s')p(s>) 



s' — s — iO 



(25) 



where gi(s') and gj(s') are vertices and p(s') is the phase space for the intermediate state. It 
is helpful to introduce the propagator matrix D^, where the non-diagonal elements D 12 = D 21 
correspond to the transitions 1 — > 2 and 2 — > 1 (see Fig. 16b). The matrix reads: 



D 



D u D 12 
D 2 \ D 22 



{Ml - s)(M 2 - s) - B l2 B 2l 
Here the following notation is used: 

M? = m*-B ii (s) i = l,2. 



Ml - s, B 12 
B 21 , Ml - s 



(26) 



(27) 



Zeros in the denominator of the propagator matrix (26) define the complex resonance masses 
after the mixing: 

n(s) = (Ml - s)(M| - s) - B 12 B 21 = . (28) 

Let us denote the complex masses of mixed states as Ma and M B . 

Consider a simple model, where the s-dependence of the function B^s) near the points 
s ~ M\ and s ~ M 2 B is assumed to be negligible. Let Ml and B 12 be constants. Then one has: 



A,B 



-(Ml + Mj) ± J -(Ml - Mf) 2 + B 12 B 21 



(29) 



In the case, when the widths of initial resonances 1 and 2 are small (hence the imaginary part 
of the transition diagram B i2 is also small), the equation (29) turns into the standard formula 
of quantum mechanics for the split of mixing levels, which become repulsive as a result of the 
mixing. If so, 



D 



cos 2 6 /(Ml -s) + sin 2 9/(M 2 B - s) - cos 9 sin 0/(M 2 A - s) + sin#cos#/(M| - s) 

cosQsm6/(M 2 A - s) + sin# cos9/(M B - s) sin 2 #/(M| - s) + cos 2 9/(M B - s) 



cos 2 9 = \ + - 



\(M 2 -M 2 ) 



^\(Ml-M 2 ) 2 + B 12 B 21 



(30) 
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Figure 16: Diagrams describing the propagation functions Du (a) and D\ 2 (b) for the interac- 
tion of two resonance states. 



The states \A > and \B > are superpositions of the initial levels, |1 > and |2 >, as follows: 



\A >= cos0|l > -sin0|2 > , \B >= sin0|l > + cos6>|2 > 



(31) 



In general, the representation of states \A > and \B > as superpositions of initial states is valid, 
when the s-dependence of functions Bij(s) can not be neglected, and their imaginary parts are 
not small. Consider the propagator matrix near s = M\: 



D = 



M|(s) - s B 12 {s) 
B 21 (s) M*(s)-s 



H'{Ml){M\-s) 



Ml{M\)-M\ B 12 (M\) 
B 2l {M\) Ml{Ml)-M\ 



(32) 

In the left-hand side of Eq. (32), only the singular (pole) terms survive. The matrix determinant 
in the right-hand side of (32) equals zero: 



[Ml{M\) - M\][Ml{Ml) - M\\ - B 12 (M 2 A )B 21 (M 2 A ) = , 



(33) 



This equality follows from Eq. (28), which fixes IT(M|) = 0. It allows us to introduce the 
complex mixing angle: 

\A >= cos0 A |l > -sin#4|2 > . (34) 

The right-hand side of Eq. (30) can be rewritten by making use of the mixing angle 9 A , as 
follows: 



D 



Mi 



cos 2 9 a 
sin 9a cos 9a sin 2 9 



— cos 9 a sin 9a 



(35) 
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where 



Na 



n'(Mi) 



[2M\- Ml - M 2 2 ], cos 2 #4 



M\ - Mf 



2.U 2 - Mf - M[ 



-, sin 2 ^ 



M 2 



Ml 



■1M\ - M( - Mf 



(36) 

We remind that in the formula (36) the functions M 2 (s), Mf (s) and 5i 2 (s) are fixed in the 
point s = M\. In the case under consideration, when the angle 9 a is a complex quantity, the 
values cos 2 #a and sin 2 6a do not determine the probability of states |1 > and |2 > in \A >; 
indeed, the values V Na cos 9a and — V Na sin 9a are the transition amplitudes \A >^ |1 > and 
| A >^ |2 >. Therefore, the corresponding probabilities are equal to | cos9a\ 2 and | sin ^a| 2 - 

In order to analyse the content of the state \B >, an analogous expansion of the propagator 
matrix should be carried out near the point s = M|. Introducing 



\B >= sin ^(l > + cos 9 B \2 > , 



(37) 



we have the following expression for D in the vicinity of the second pole s = M^: 



D 



'M% M| - s 



sm 2 9 B cos 9 B sin 9 B 

sin 9b cos 9 b cos 2 9 b 



(38) 



where 



N. 



n'(M 2 ) 



2Mf - Ml - Ml , cos 2 9 B 



M 2 B - Ml 



2M| - Ml - Ml 



, sin 2 9 B 



M 2 B - Ml 



2Ml - Ml - Ml ' 
(39) 

In Eqs. (38), (39) the functions M 2 (s), Mf(s) and B 12 (s) are fixed in the point s = M 2 B . 

If B\2 depends on s weakly and one can neglect this dependence, the angles 9a and 9 B 
coincide. In general, however, they are different. So the formulae for the propagator matrix 
differ from the standard approach of quantum mechanics by this very point. 

Another distinction is related to the type of the level shift afforded by mixing, namely, in 
quantum mechanics the levels "repulse" each other from the mean value l/2(Ei + E 2 ) (see 
also Eq. (29)). Generally speaking, the equation (28) can cause both a "repulsion" of masses 
squared from the mean value, 1/2(M 2 + Mf), and an "attraction". 

Let us remind now, how to write the amplitudes in the one-channel and multi-channel cases. 
The scattering amplitude in the one-channel case is defined by the following expression: 



M s ) = 5 , i( s )Aj(s)fi'j(^ 



(40) 



u 



In the multi-channel case, Bij(s) is a sum of loop diagrams: Bij(s) = J2 a Bif(s), where B\ a ^ 

is a loop diagram in the channel a with vertex functions g^ a \ 
The partial scattering amplitude in the channel a — > b equals 



and a phase space factor p a . 



Aa^(s) = g < i a) (s)D l ,(s)g) >(s). 



(41) 
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5.2.2 Construction of propagator matrix in a general case (N resonances) 



Consider the construction of the propagator matrix D for an arbitrary number (N) of reso- 
nances. The matrix elements, D^, describe the transition from the initial state % (with the bare 
propagator {mf — s)~ v ) to the state j. They obey the system of linear equations as follows: 



Dij = D ik B kj (s)(rr£ - s) 1 + <^(mj - s) 



-i 



(42) 



where B^s) is the loop diagram for the transition % — > j and Sij is the Kronecker symbol. 
Let us introduce the diagonal propagator matrix d for initial states : 

d = diag {{rn\ — s) _1 , (m\ — s) _1 , (m\ — s)" 1 • • •) ■ 

Then the system of linear equations (39) can be rewritten in the matrix form as 

D = DBd + d . 

One obtains 

T 

D = 



(d-i - B) 

The matrix d~ l is diagonal, hence D^ 1 = (d^ 1 — B) is of the form 

Ml-s -B 12 (s) -B 13 (s) ■■■ 



(43) 

(44) 
(45) 



b- 1 



-B 21 (s) Ml-s -B 23 (s) 
-B 31 (s) -B 32 (s) Ml-s 



(46) 



where Mf is defined by Eq. (27). Inverting this matrix, we obtain a full set of elements Dij(s): 

(47) 



lj{s> nw( s ) 



Here n^(s) is the determinant of the matrix D 1 , and Ir^ '(s) is a matrix supplement to 
the element [-D -1 ]^, i.e. the matrix D^ 1 with an excluded j-th line and i-th column. 

The zeros of n^(s) define the poles of the propagator matrix which correspond to physical 
resonances formed by the mixing. We denote the complex resonance masses as: 



s = M*, Ml, M 2 C ,... 



(48) 



Near the point s = M\, one can leave in the propagator matrix the leading pole term only. 
This means that the free term in Eq. (44) can be neglected, so we get a system of homogeneous 
equations: 

D lk (s)(d- 1 -B) k .=0. (49) 
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The solution of this system is defined up to the normalisation factor, and it does not depend on 
the initial index i. If so, the elements of the propagator matrix may be written in a factorised 
form as follows: 



N A 



of, ai«2) cxia^, 
a 2 ai, a%, ot 2 a 3 , 
0:301, 0:30:2, CK3, 



where Na is the normalisation factor chosen to satisfy the condition 



9 9 



a, 



a 



N 



1. 



(50) 



(51) 



The constants Oj are the normalised amplitudes for the resonance A — > state i transitions. The 
probability to find the state i in the physical resonance A is equal to: 



Wi = Oti 



(52) 



Analogous representations of the propagator matrix can be given also in the vicinity of other 
poles: 



D\f\ S ~ Ml) = /V BJ/2 



(53) 



The coupling constants satisfy normalisation conditions similar to that of Eq. (48): 



7l 2 + 7 2 2 + . . . + 7^ = 1 



(54) 



In the general case, however, there is no completeness condition for the inverse expansion: 

«• +A 2 + 7, 2 + ...^l. (55) 

For two resonances this means that cos 2 0^ + sin 2 9g ^ 1. Still, let us remind that the equality 
in the inverse expansion, which is relevant for the completeness condition, appears in models 
where the s-dependence of the loop diagrams is neglected. 



5.2.3 Full resonance overlapping: the accumulation of widths of neighbouring 
resonances by one of them 

Let us consider two examples which describe the idealised situation of a full overlapping of two 
or three resonances. In these examples, the effect of accumulation of widths of neighbouring 
resonances by one of them can be seen in its original untouched form. 

a) Full overlapping of two resonances. 

For the sake of simplicity, let Bij be a weak s-dependent function; Eq. (29) may be used. 
We define: 

Ml = M 2 - iM R T 1 , Ml = M 2 - iM*r 2 , (56) 

and put 

Re Bl2 (Mi)=P / ^fy^ G. (57) 



2 



71 s' - M; 
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It is possible that Rei?i2(M^) equals zero at positive g\ and g 2 , if the contribution from the 
integration region s' < M R cancels the contribution from the s' > M R region. In this case 

B 12 (M 2 R ) -> ig 1 (M 2 R )g 2 (M 2 R )p(M 2 R ) = iM R J^T~ 2 . (58) 

Substituting Eqs. (53)-(55) into Eq. (45), one has: 

M 2 A ^M 2 R -tM R (T 1 + T 2 ) M 2 B ~^M 2 R . (59) 

Therefore, after mixing, one of the states accumulates the widths of primary resonances, Ta — > 
Ti + T 2 , and another state becomes a quasi-stable particle, with T B — > 0. 

b) Full overlapping of three resonances. 

Consider the equation 

n^(s) = o (60) 

in the same approximation as in the above example. Correspondingly, we put: 

ReB ab (M R ) ->0,(a^6); M 2 = M 2 R - s - iM R Y t = x - t% . (61) 

A new variable, x = M R — s, is used, and we denote M R Yi = 7,. Taking into account BijBji = 
— 7i7j and Bi 2 B 23 B 3 i = —2717273, we can rewrite the equation (60) as follows: 

x 3 + ar 2 (i7i + 172 + 173) = . (62) 

Therefore, at full resonance overlapping, one obtains: 

M\ -> M 2 R - iM R (T 1 + T 2 + r 3 ) , M 2 B ^M 2 R , M 2 -^M 2 R . (63) 

Thus, the resonance A has accumulated the widths of three primary resonances, and the states 
B and C became quasi-stable and degenerate. 



6 The qq-gg content of /2-niesons, observed in the reac- 
tions pp — ► 7r°7r°, 7777, r]r}' 

We determine here the qq—gg content of / 2 -mesons, observed in the reactions pp — > 7r°7r°, 7777, 7777' 
[2] . This determination is based on experimentally observed relations (3) and the rules of quark 
combinatorics taken into account in the leading terms of the 1/iV-expansion. 

For the f 2 — > 7T 7r°, 7777, 7777' transitions, when the qq-meson is a mixture of quarkonium and 
gluonium components, the decay vertices read in the leading terms of the 1/iV-expansion (see 
Tables 3 and 4) as follows: 

qq— meson COS G 



grt-meson = g f 2 Q C _^V + ^2 Q^ sin ^ + G ( cos ^ Q + \ sin 2 ©) ) 

V v2 / v2 + A 



qq— meson ■ r\ r\ 

g™, = sin 6 cos 6 



f COS if rr . \ G 



(64) 
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The terms proportional to g stand for the qq — > two mesons transitions (g = g L cosa), while 
the terms with G represent the gluonium — > two mesons transition (G = G L sin a) . Conse- 
quently, G 2 and g 2 are proportional to the probabilities for finding gluonium (W = sin 2 a) and 
quarkonium (1 — W) components in the considered / 2 -meson. Let us remind that the mixing 
angle stands for the nn and ss components in the n and rj' mesons; we neglect the possible 
admixture of a gluonium component to n and n' (according to [47] , the gluonium admixture to 
n is less than 5%, to rj' — less than 20%). For the mixing angle 6 we take 6 = 37°. 



6.1 The analysis of the quarkonium-gluonium contents of 
the / 2 (1920), / 2 (2020), / 2 (2240), / 2 (2300) 

Making use of the data (3), the relations (64) allow us to to find if as a function of the ratio G/g 
of the coupling constants. The result for the resonances / 2 (1920), / 2 (2020), / 2 (2240), / 2 (2300) 
is shown in Fig. 17. Solid curves enclose the values of g m j 'g^o^o for A = 0.6 (this is the 7/77-zone 
in the (G/g,<p) plane) and dashed curves enclose g m i / g^o for A = 0.6 (the nn'-zone). The 
values of G/g and if, lying in both zones, describe the experimental data (3): these regions are 
shadowed in Fig. 17. 

The correlation curves in Fig. 17 enable us to give a qualitative estimate for the change of 
the angle ip (i.e. the relation of the nn and ss components in the / 2 meson) depending on the 
value of the gluonium admixture. The values g 2 and G 2 are proportional to the probabilities 
of having quarkonium and gluonium components in the / 2 meson, g 2 = (g L ) 2 (l — W) and 
G 2 = (G L ) 2 W . Here W is the probability of a gluonium admixture in the considered gg-meson; 
g L and G L are universal constants, see Tables 3 and 4. Since G L /g L ~ \/y/W c and W ~ l/N c , 
we can give a rough estimate: 

g~ W ^ — . (65) 
(f N c (l-W) 10 V ' 

Let us remind that the numerical calculations of the diagrams indicate that 1/N C leads to a 
smallness of the order of 1/10 - this is taken into account in (65). Assuming that the gluonium 
components are less than 20% (W < 0.2) in each of the qq resonances / 2 (1920), / 2 (2020), 
/ 2 (2240), / 2 (2300), we put roughly W ~ 10 G 2 /g 2 , and obtain for the angles (p the following 
intervals: 



W gluoniurn [f 2 (1920)] < 20% 
W gluomum [f 2 (2020)] < 20% 
W gluomum [f 2 (22A0)] < 20% 
H T gluonium [/ 2 (2300)] < 20% 



-0.8° < y>[/ 2 (1920)] < 3.6° , 

-7.5° < y>[/ 2 (2020)] < 13.2° , 

-8.3° < v?[/ 2 (2240)] < 17.3° , 

-25.6° < v?[/ 2 (2300)] < 9.3° (66) 



6.2 The nn-ss content of the qrg-mesons 

Let us summarise what we know about the status of the (I = 0, J PC = 2 ++ ) gg-mesons. 
Estimating the nn-ss content of the / 2 -mesons, we ignore the gg admixture (remembering that 
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it is of the order of sin 2 a ~ l/N c ). 

1. The resonances / 2 (1270) and /2(1525) are well-known partners of the basic nonet with 
n = 1 and a dominant P-component, l 3 P 2 gg. Their flavour content, obtained from the 
reaction 77 — > K$Ks, is 

/ 2 (1270) = cos(/9 n= inn + sin(/9 n=1 ss, 
/ 2 (1525) = — sin ip n= inn + cos ip n= iss, 

<p n=1 = -1±3°. (67) 

2. The resonances / 2 (1560) and / 2 (1750) are partners in a nonet with n = 2 and a dominant 
P-component, 2 3 P 2 qq. Their flavour content, obtained from the reaction 77 — > K S K S) is 

/ 2 (1560) = cosip n=2 nn + siny2 n=2 ss, 
/ 2 (1750) = - sin (p n=2 nn + cos y?n=2ss, 

^„=i = -10lf °. (68) 

3. The resonances / 2 (1920) and / 2 (2120) [4] (in [1] they are denoted as / 2 (1910) and / 2 (2010)) 
are partners in a nonet with n = 3 and with a dominant P-component, 3 3 P 2 gg. Ignor- 
ing the contribution of the glueball component, their flavour content, obtained from the 
reactions pp — > 7r°7r°, r/r], rjr)', is 

/ 2 (1920) = cos y? n=3 nn + sin ip n=3 ss, 
/ 2 (2120) = - sin (p n=3 nn + cos y? n=3 ss, 

ip n=3 = 0±5°. (69) 

4. The next, predominantly 3 P 2 states with n = 4 are / 2 (2240) and / 2 (2410) [4]. (By mistake, 
in [1] the resonance / 2 (2240) [2] is listed as / 2 (2300), while / 2 (2410) [4] is denoted as 
/ 2 (2340)). Their flavour content at W = is determined as 

/ 2 (2240) = cos(p n= 4nn + sin <^ n=4 ss, 
/ 2 (2410) = — sin (p n=i nn + cos (fn^ss, 

V? n=4 = 5 ±11°. (70) 

5. / 2 (2020) and / 2 (2340) [4] belong to the basic P-wave nonet (n = 1) (in [1] the / 2 (2020) 
[2] is denoted as / 2 (2000) and is put in the section "Other light mesons", while / 2 (2340) 
[4] is denoted as / 2 (2300)). The flavour content of the 1 3 P 2 mesons is 
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/ 2 (2020) = cos(p n ( F ) =1 nn + sin (p n ( F ) =1 ss, 
/ 2 (2340) = - sin (p nlyF)=1 nn + cos (p n (F)=iSS, 
ip n (F)=i = 5 ± 8°. 



(71) 



6. The resonance / 2 (2300) [2] has a dominant F-wave quark- ant iquark component; its flavour 
content for W = is defined as 



A partner of / 2 (2300) in the 2 3 F 2 nonet has to be a / 2 -resonance with a mass M ~ 
2570 MeV. 



7 Conclusion 

The broad / 2 (2000) state is the descendant of the lowest tensor glueball. This statement is 
favoured by estimates of parameters of the Pomeron trajectory (e.g., see [5], Chapter 5.4, 
and references therein), according to which M 2 ++ g i ue bau — 1.7 — 2.5 GeV. Lattice calculations 
result in a similar value, namely, 2.2-2.4 GeV [21]. The corresponding coupling constants 
/ 2 (2000) — > 7T°iT ,riri,riri' satisfy the relations for the glueball, eq.(22), with A ~ 0.5 — 0.7. The 
admixture of the quarkonium component (qq) g i U ebaU in ^2(2000) cannot be determined by the 
ratios of the coupling constants between the hadronic channels; to define it, / 2 (2000) has to be 
observed in 77-collisions. The value of (qq) g i U ebaii in /2(2000) may be rather large: the rules 
of 1 /iV-expansion give a value of the order of Nf/N c . It is, probably, just the largeness of 
the quark-antiquark component in / 2 (2000) which results in its suppressed production in the 
radiative J ftp decays (see discussion in [8]). 

We have now two observed glueballs, a scalar / (1200 — 1600) [14, 15] (see also [5, 12]) 
and a tensor one, / 2 (2000). It is illustrative to present the situation with ++ ,2 ++ glueballs 
on the (J, M 2 )-plane, we demonstrate this in Fig. 18. According to various estimates, the 
leading Pomeron trajectory has an intercept at a(0) ~ 1.10 — 1.30 (see, for example, [18, 19, 
20]). Assuming that the Pomeron trajectory has a linear behaviour, which does not contradict 
experimental data, a P (M 2 ) = a P (0) + a' P (0)M 2 , we have for the slope a' P (0) = 0.20 ± 0.05. 
The scalar glueball / (1200 — 1600) is located on the daughter trajectory which predicts the 
second tensor glueball at M ~ 3.45 GeV. If the Pomeron trajectories in (n, M 2 ) plane are linear, 
similar to the gg-trajectories, then the next scalar glueball (radial exitation of gg gluonium) 
should be at M ~ 3.2Ge\/. 

Observed glueball states have transformed into broad resonances owing to the accumulation 
of widths of their neighbours. The existence of a low-lying pseudoscalar glueball is also expected. 
It is natural to assume that it has also transformed into a broad resonance. Consequently, the 
question is, where to look for this broad _+ state. There are two regions in which we can 



/ 2 (2300) = cos <f n (F)=2nn + sin (p n ( F )=2SS, (p n (F)=2 = -8° ± 12°. 



(72) 
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suspect the existence of a pseudoscalar glueball: in the region of 1700 MeV or much higher, 
at ~ 2300 MeV, see the discussion in [8] (Section 10.5). In [50] it is suggested that the lowest 
scalar and pseudoscalar glueballs must have roughly equal masses. If so, a 0~ + glueball has to 
occur in the 1700 MeV region. 

The authors are grateful to D.V. Bugg, L.D. Faddeev and S.S. Gershtein for stimulating 
discussions. The paper was supported by the grant No. 04-02-17091 of the RFFI. 



40 



References 

[1] S. Eidelman et al (PDG), Phys. Lett. B 592, 1 (2004). 

[2] A.V. Anisovich et al, Phys. Lett. B 491, 47 (2000). 

[3] D. Barberis et al. (WA 102 Collab.), Phys. Lett. B 471, 440 (2000). 

[4] R.S. Longacre and S.J. Lindenbaum, Report BNL- 72371-2004. 

[5] V.V. Anisovich, M.N. Kobrinsky, J. Nyiri, Yu.M. Shabelski, "Quark model and high energy 
collisions" , second edition, World Scientific, 2004. 

[6] V.A. Schegelsky, A.V. Sarantsev, V.A. Nikonov, " Phenomenological investigation of the 
K S K S final state in two-photon collisions and nonet classification of tensor resonances", 
L3 Note 3001, October 27, 2004. 

[7] A. Etkin et al, Phys. Lett. B 165, 217 (1985); B 201, 568 (1988). 

[8] D.V. Bugg, Phys. Rep., 397, 257 (2004). 

[9] V.V. Anisovich, Pis'ma v ZhETF, 80, 845 (2004) [JETP Letters, 80, 715 (2004)]. 

[10] V.A. Schegelsky and A.V. Sarantsev, Resonanses in 77 — > 3n reaction, Talk given at 
XXXIX PNPI Winter School (2005). 

[11] A.V. Anisovich, V.V. Anisovich and A.V. Sarantsev, Phys. Rev. D 62, 051502 (2000). 

[12] V.V. Anisovich, UFN, 174, 49 (2004) [Physics-Uspekhi, 47, 45 (2004)]. 

[13] V.V. Anisovich, D.V. Bugg and A.V. Sarantsev, Phys. Rev. D 58, 111503 (1998). 

[14] V.V. Anisovich, Yu.D. Prokoshkin and A.V. Sarantsev, Phys. Lett. B 389, 388 (1996); Z. 
Phys. A 357, 123 (1997). 

[15] A.V. Anisovich, V.V. Anisovich and A.V. Sarantsev, Phys. Lett. B 395, 123 (1997); Z. 
Phys. A 359, 173 (1997). 

[16] V.V. Anisovich, M.A. Matveev, J. Nyiri and A.V. Sarantsev, Yad. Fiz., in press hep- 
ph/0501003. 

[17] V.V. Anisovich and A.V. Sarantsev, Pis'ma v ZhETF, 81,531 (2005), [JETP Letters, 81, 
417 (2005)]. 

[18] A.B. Kaidalov and K.A. Ter-Martirosyan, Sov. J. Nucl. Phys. 39, 979 (1984). 

[19] P.V. Landshoff, "Soft hadron reactions", in QCD: 20 Years Later, eds. P.M. Zerwas and 
H.A. Kastrup, (World Scientific, Singapore, 1993). 

[20] L.G. Dakhno and V.A. Nikonov, Eur. Phys. J. A5, 209 (1999). 

41 



[21] G.S. Bali, K. Schilling, A. Hulsebos et al. (UK QCD Collab.), Phys. Lett. B 309, 378 
(1993); 

C.J. Morningstar, M.J. Peardon, Phys. Rev. D 60, 034509 
M. Loan, X-Q. Luo and Z-H.Luo, hep-lat/0503038. (1999). 

[22] U. Loering, B.C. Metsch and H.R. Petry, Eur. Phys. J. A10, 447 (2001). 

[23] E. Klempt, Eur.Phys.J. C28 (2003), hep-ex/0206012. 

[24] D.V. Bugg and A.V. Sarantsev, B.S. Zou, Nucl. Phys. B471, 59 (1996). 

[25] G. Grayer et al., Nucl. Phys. B75, 189 (1974). 

[26] V.V. Anisovich and A.V. Sarantsev, Eur. Phys. J. A 16, 229 (2003). 

[27] D.V. Amelin et al. (VES Collab.), Z. Phys. C70, 71 (1996). 

[28] V.V. Anisovich, D.I. Melikhov, V.A. Nikonov, Phys. Rev. D 55, 2918 (1997). 

[29] A.V. Anisovich, V.V. Anisovich and V.A. Nikonov, Eur. Phys. J. A12, 103 (2001); 

A.V. Anisovich, V.V. Anisovich, M.A. Matveev and V.A. Nikonov, Yad. Fiz. 66, 946 (2003) 
[Phys. Atom. Nucl. 66, 914 (2003)]. 

[30] A.V. Anisovich, V.V. Anisovich, L.G. Dakhno, V.A. Nikonov and A.V. Sarantsev, Yad. 
Fiz., in press, hep-ph/0406320. 

[31] J.Z. Bai et al., Phys. Rev. D68, 052003 (2003); Phys. Lett. B 472, 207 (2000). 

[32] V.V. Anisovich, UFN, 168, 481 (1998) [Physics-Uspekhi, 41, 419 (1998)]. 

[33] G. 't Hooft, Nucl. Phys. B 72, 461 (1974). 

[34] G. Veneziano, Nucl. Phys. B 117, 519 (1976). 

[35] F.E. Close and A. Kirk, Eur. Phys. J. C21, 531 (2001). 

[36] W. Lee and D. Weingarten, Phys. Rev. D60, 034509 (1999). 

[37] I.S. Shapiro, Nucl. Phys. A 122 645 (1968). 

[38] I.Yu. Kobzarev, N.N. Nikolaev, L.B. Okun, Yad. Fiz. 10, 864 (1969); [Sov. J. Nucl. Phys. 
10, 499 (I960)]. 

[39] L. Stodolsky, Phys. Rev. D 1, 2683 (1970). 

[40] A.V. Anisovich, V.A. Nikonov, A.V. Sarantsev, V.V. Sarantsev, in "PNPI XXX, 
Scientific Highlight, Theoretical Physics Division, Gatchina (2001), p. 58.(1999). 

[41] S.M. Flatte, Phys. Lett. B63, 224 (1976). 



42 



[42] E. Eisenhandler et al, Nucl. Phys. B98, 109 (1975); 
A. Hasan et al., Nucl. Phys. B378, 3 (1992). 

[43] A.V. Anisovich, hep-ph/0104005. 

[44] V.V. Anisovich, M.G. Huber, M.N. Kobrinsky and B.Ch. Metsch, Phys. Rev. D 42, 3045 
(1990); 

V.V. Anisovich, V.A. Nikonov and J. Nyiri, Yad. Fiz. 64, 1 (2001). 
[45] K. Peters, E. Klempt, Phys. Lett. B 352, 467 (1995). 

[46] S.S. Gershtein, A.K. Likhoded, Yu.D. Prokoshkin, Z. Phys. C 24, 305 (1984); 

C. Amsler, F.E. Close, Phys. Rev. D 53, 295 (1996); Phys. Lett. B 353, 385 (1995); 
V.V. Anisovich, Phys. Lett. B 364, 195 (1995). 

[47] V.V. Anisovich, D.V. Bugg, D.I. Melikhov, V.A. Nikonov, Phys. Lett. B 404, 166 (1997). 

[48] A.A. Carter et al, Phys. Lett. B 67, 117 (1977). 

[49] A.V. Anisovich, V.V. Anisovich, V.N. Markov, M.A. Matveev and A.V. Sarantsev, J. Phys. 
G: Nucl. Part. Phys. 28, 15 (2002). 

[50] L. Faddeev, A.J. Niemi and U. Wiedner, hep-ph/0308240 (2003). 

[51] A.V. Anisovich et al, Phys. Lett. B 452, 173 (1999); Phys. Lett. B 452, 187 (1999); Phys. 
Lett. B 517, 261 (2001); 



43 




Figure 17: Correlation curves g rm ((p,G/g)/g n o n o(ip,G/g) and g m * (<p, G / g) / g n o n o (ip, G / g) drawn 
according to (64) at A = 0.6 for / 2 (1920), / 2 (2020), / 2 (2240), / 2 (2300) [2,17]. Solid and dashed 
curves enclose the values g^i^, G / ' g) / 'g n o n o(ip, G/g) and g m '(<f, G/g)/ g n o n o(ip, G/g) which obey 
(3) (the zones 7777 and 777/ in the {G/g, ip) plane). The values of G/g and ip, lying in both zones 
describe the experimental data (3): these are the shadowed regions. 
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Figure 18: Glueball states on the Pomeron trajectories. 
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